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Abstract. Let <f> be an endomorphism of the projective line of degree at 
least 2, defined over a noetherian commutative ring R with unity. We show 
that the automorphism group of <j> is a finite group scheme, and we construct 
algorithms to compute it when R is a finite field or a number field. We also give 
an algorithm for determining when two such endomorphisms are conjugate. We 
have implemented these algorithms in Sage when R is a finite field or the field 
of rational numbers. 



1. Introduction 

Let F be a field, and let <f> = f/g e F(z) be a rational function such that 
gcd(/, g) = 1 and d = deg(0) := max{deg(/), deg(g)} > I. When viewed as an 
endomorphism of the projective line <fi : P^ — > P^, a dynamical theory of <f> arises 
from iteration. That is, for x € P^, we may consider its orbit 

x n- cj)(x) i — ^ (j) 2 {x) n- (j) z (x) i-> • • • 

(Here we write (f) 1 = <fi and (j> n = <j> o (ff 1 ^ 1 for each n > 1.) The case F = 
C — dynamics of self-maps of the Riemann sphere — has a fascinating history 
dating back as far as Newton; e.g., see [TJ [TT]. When F is a finite field, these 
dynamical systems behave (conjecturally) like random maps, which has applications 
to factoring integers [131 [5] . If F is a number field, then we have the younger theory 
of arithmetic dynamics [15) . The case where F is a non- Archimedean field is younger 
still and draws much inspiration from the complex case [U [10] . 

In the present paper, we study algorithmic aspects of a fundamental question: 

When do two rational functions <f>, ip € F(z) exhibit the same 
dynamical behavior? 

For topological reasons, they must have the same degree. Given any fractional linear 
transformation f(z) = ^f^f , viewed as an clement of PGL2(F), the dynamical 
behavior of <j> is equivalent to that of <fy := f o <fi o indeed, / maps the </>-orbit 
of a point x G P^ to the <^-orbit of f(x). If there exists / G PGL2(F) such that 
ip = f o <f> o then we say that <j> and ifi are conjugate (over F). Conjugation 
defines an action of PGL2(F) on the parameter space of rational functions of fixed 
degree d > 1, denoted Katd(F). The stabilizer of a rational function <j> for this 
action is called the automorphism group of </>, and is denoted by Aut^(F). This 
group is always finite, and it is trivial for most rational functions. 
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We thus have two problems which happen to be computationally quite similar: 
(1) determine whether two given rational functions are conjugate, and (2) determine 
whether a given rational function has nontrivial automorphism group. We focus 
mainly on algorithms designed to compute automorphism groups, and in the final 
section we sketch the modifications needed to address the first problem. 

We design three algorithms for computing automorphism groups of rational func- 
tions, each applying to a slightly different setting. Let F be a field and <fi : Pp —> Pp 
be a morphism of degree at least 2, as above. Our first algorithm computes the 
absolute automorphism group Aut^-F) and a field of definition E/F; that is, 
Aut^(E) = Aut^(F). This algorithm requires constructing the splitting field of 
a polynomial with degree 0(d), so this is not very practical over number fields un- 
less d is small. Over finite fields it is much more efficient because the average size 
of the splitting field of a polynomial of degree 0(d) is significantly smaller. 

For s e PGL 2 (F), write Fix(s) for the set of fixed points of s in P 1 ^). If 
s is an automorphism of cf>, the action of <j) on Fix(s) is highly restricted, both 
geometrically and arithmetically. Our second algorithm takes advantage of this 
fact to compute Aut^(F) when F is any field of characteristic or when F is a 
finite field. If F has characteristic p > and is not finite, then the algorithm only 
detects the elements of Aut^F) whose order is prime to the characteristic. This 
algorithm requires finding linear and quadratic factors of a polynomial of degree 
d 2 + 1. With the present implementations of root finding and polynomial factoring 
over number fields available in Sage and Magma, this is infeasible when d is large. 
However, the algorithm is quite efficient when F is a finite field and for number 
fields when d is reasonable, say d < 15. 

Our third algorithm computes Aut,p(F) when F is a number field. For rational 
maps of degree at least 15 it is significantly faster than the second method. The 
main idea here is to first reduce the coefficients modulo v for several finite places 
v of good reduction and compute the automorphism group over the residue field 
(using our second algorithm) . We then use the Chinese remainder theorem to piece 
together the various automorphisms modulo v to arrive at the set of automorphisms 
over the original number field F. 

We first prove a general result about the algebro- geometric structure of Aut^. 
As we will work with automorphism groups over global fields and the reductions at 
several primes, this type of result is necessary to ensure we are on solid footing. For 
notation, let R be a noetherian commutative ring with unity, and let i?-Alg and 
Grp denote the categories of commutative i?-algebras and groups, respectively. For 
any i?-algebra S, we identify PGI^jS 1 ) with Aut(Pg), the group of automorphisms 
of P 1 defined over S. We make the following definition: 

Definition. Let <f> : Pp —> Pp be a nonconstant morphism. The automorphism 

group of <fi is the i?-group scheme Aut^ represented by the functor Aut^ : i?-Alg — s- 
Grp defined by 

Aut^(S) = {/ G Aut(P^) : <!> = := / o <f> o J" 1 }. 

Theorem 1.1. Let R be a commutative ring and let <f> : Pp —> Pp be a noncon- 
stant endomorphism. Then the functor Aut ^ is represented by a closed R-subgroup 
scheme Aut^ C PGL2. If moreover deg(<^) > 2, then Aut^ is finite over Speci?. 

Remark 1.2. The group scheme Aut^ need not be flat over Speci?. For example, 
if 4>(z) = z 1 as an endomorphism of P^, then one can check that Aut^(Q) = 
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{z,l/z} = Aut (F p ) for p > 2. But Aut (F 2 ) ^ PGL 2 (F 2 ), which has order 6. 
Since the order of a finite flat group scheme is locally constant, we conclude that 
Aut^ is not flat over Spec Z. 

Remark f.3. If deg(4>) = 1, then Aut0 is not a finite group scheme in general. 
Consider the examples 4>(z) = z and ip(z) = z + 1, for which Aut^ = PGL 2 and 
Aut^ = G a , respectively. 

The next result will have two uses in this paper. First, it will allow us to deduce 
that Aut^, is proper when <f> has degree at least 2. Second, it will provide the main 
tool for relating the automorphism group of an endomorphism over a number field 
to the automorphism group over a finite field. For the statement, if k is a non- 
Archimedean field (not necessarily complete) with valuation ring o, we say that 
an endomorphism <j> : — > has good reduction if there exists a morphism 
$ : Pj — >• Pj with generic fiber <f>. 

Reduction Lemma. Let k be a non- Archimedean field with valuation ring o and 
residue field F, and let <j) G k(z) be a rational function of degree at least 2 (which 
is equivalent to a morphism P£ — > Pi). Suppose that <f> has good reduction. Then 
every element of Aut^(fc) has good reduction, and the canonical reduction o — > F 
induces a homomorphism red : Aut^(fc) — > Aut^(F). IfW has characteristic p > 
(resp. characteristic zero), then the kernel of reduction is a p-group (resp. trivial). 

If K is a number field and v is a finite place of K, we write K Vl Z v , and F„ for 
the completion of K at v, the valuation ring of K v , and the residue field of K V1 
respectively. If G K(z) is a rational function, we say that it has good reduction 
at v if all of its coefficients arc integral at v, and it extends to a morphism over 
SpecZ„. (Equivalcntly, <f> has good reduction if one can reduce its coefficients 
modulo v, and the resulting morphism of Pp has the same degree as cf>. ) 

Proposition 1.4. Let K be a number field and let <j) <= K(z) a rational function of 
degree d > 2. Define Sq to be the set of rational primes given by 



and let S be the (finite) set of places of K of bad reduction for <f> along with the 
places that divide a prime in Sq. Then red.„ : Ant <p(K) — > Aut^(F„) is a well-defined 
injective homomorphism for all places v outside S. 

Remark 1.5. In practice, Proposition 11.41 allows one to determine the group struc- 
ture of Aut,p(K) very quickly by computing Aut^(F„) for a few places v S. This 
is analogous to the way one typically computes the torsion subgroup of an elliptic 
curve over a number field; see [EJ VII. 3]. If one wishes to compute the elements of 
Aut^(A') rather than just the group structure, then more work is required. 

Outline. Sectionals occupied by the proof of the Reduction Lemma and its corol- 
lary. The argument uses techniques from the ramification theory of endomorphisms 
of the Berkovich projective line; the reader may safely skip the proof and use the 
Reduction Lemma as a black box if necessary. The proof of Theorem 11.11 is given 
in Section |3l Section [4] is occupied by the algorithms for computing automorphism 
groups. Section [5] contains a few examples of our implementation of the algorithms. 
All of the computations were carried out using Sage |17j . Our code is included with 
the arXiv distribution of this article. Finally, Section |6] contains a brief discussion 
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of the scheme Conj^, ^ of elements of PGL2 that conjugate <f> to ip, and we also 
sketch the main ideas of an algorithm for determining if two endomorphisms </> and 
ip are conjugate, i.e., if the scheme Conj^ ^ has any points. 
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2. Proof of the Reduction Lemma 

For background on the Berkovich projective line and dynamics, see [3]. For a 
more concise summary of the necessary ideas, we direct the reader to [BJ. 

Let Cfc be the completion of an algebraic closure of the completion of fc, and 
let P 1 be the Berkovich analytification of the projective line P^ fc . The morphism 
(j> extends functorially to P 1 . We use two key facts due to Rivera-Letelier [14j 
Thm. 4]: (1) <f> has good reduction if and only if the Gauss point ( e P 1 is totally 
invariant, and (2) a rational function has at most one totally invariant point in 
P 1 \¥\C k ). 

For / £ Aut ( p(k), we have 

r\o = /- 1 (r 1 (0) = w>° f)-\o = (/ ° <t>r\o = r\r\o)- 

Hence / -1 (C) is a totally invariant type II point for 0, so that /(C) = C- Equiva- 
lent^, / has good reduction. Thus the reduction map red : Aut^(fc) — > Aut^(F) is 
well-defined, and it is evidently a homomorphism. 

Now we compute the kernel of reduction. Suppose red(/) is trivial. Without loss 
of generality, we may replace k with a finite extension in order to assume that / has 
a fc-rational fixed point. Moreover, we may conjugate / by an element of PGL2(o) 
in order to assume that /(oo) = 00. Now f(z) = az + (3. If / has order m > 1, 
then the equation f m (z) = f(z) shows that a is an m-th root of unity. But red(/) 
is trivial, so we have a = 1. If k has residue characteristic zero, then we conclude 
that a = 1 and /3 = 0. Otherwise, we find that a is a p-power root of unity in fc, 
and hence / has p-power order in Aut^(fc). The proof of the Reduction Lemma is 
complete. 

Remark 2.1 . A different proof of the first part of the Reduction Lemma can be given 
using the maximum modulus principle in non- Archimedean analysis |12| Lem. 6]. 

Proposition 2.2. Let F be a field, and let n > 2 be an integer. Suppose that 
<fi : Vp — > Vp is a morphism of degree d > 2 such that Aut ,p(F) contains an element 
of order n. Then n divides d{d 2 — 1). 

Proof. We may assume without loss of generality that F is algebraically closed. Let 
s G Aut^(F) have order n. We conjugate one of the fixed points of s to 00, so that 
s = ( a i)- (Note that replacing s with usu^ 1 has the effect of replacing <p with 
wo if) o vT 1 .) The proof divides into two cases, depending on whether s has one or 
two fixed points. 
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If s has only one fixed point, then necessarily n = char(_F) is prime and a = 1. 
(See, for example, [5j Lem. 3.1].) Replace s with 1 1 ^ s ( ^ 1 ) in order to assume 

that (8 = 1. It follows that <p(z + 1) — 1 = <j>{z), or equivalently, that the function 
4>{z) — z is invariant under the map z i— > z + 1. Hence there exists a rational function 
^(z) € -FXz) such that — z = — z). We conclude that deg(</>) = n ■ deg(ip) 
or n ■ deg(ip) + 1. 

Now suppose that s has two distinct fixed points: oo and (3/(1 — a). We may 
conjugate the second fixed point to in order to assume that = 0. Note that this 
implies that a <= F x has multiplicative order n. To say that s is an automorphism of 
4> is equivalent to saying that 4>(z) / z is invariant under the map z <— > az. Hence there 
is a rational function ip £ F(z) such that <\>(z)jz = tp(z n ). So dcg(0) = n ■ deg(tp) 
or n ■ deg(?/>) ±1- □ 

Proof of Proposition \l-4\ By the Reduction Lemma, it suffices to prove that if v ^ 
S, then Aut,f,(K) has no element of order p, where v \ p. Suppose otherwise. 

The group PGL2(A') contains an element of order p if and only if ( p + Cp 1 € K 
for some primitive p-th root of unity Cp [J- Note that [Q(Cp + Cp 1 ) '■ Q] = \{p— 1) 
for p > 2, so that \ [K : Q]. If Aut^K) contains an element of order p, then p 
divides d(d 2 — 1) by Proposition 12.21 Hence p £ So, and so v € S. □ 



3. Proof of Theorem 11.11 

Fix a commutative ring R. Over R, PGL2 may be embedded as an affine subvari- 
ety of P'jj = Proj i?[a, b, c, d]; indeed, it is the complement of the quadric ad— be = 0. 
Let <f) : P^ — > P}j be a nonconstant endomorphism. We may define Aut^ as a sub- 
group scheme of PGL2 as follows. After fixing coordinates of Pjj, the morphism 
<j) can be given by a pair of homogeneous polynomials $ = ($o(A, F), $i{X, F)) 
of degree D = deg(<f>) with coefficients in R such that the homogeneous resultant 
Res(<i>o, $1) is a unit in R. The pair $01^1 is unique up to multiplication by a 
common unit in R. Similarly, for any i?-algebra S, an element / 6 PGL2(5) may 
be given by a pair F = (aX + bY, cX + dY) with a, b,c,d £ S and ad — be € S x . 
Note that Z^ 1 is represented by the pair F^ 1 := (dA — 6F, — cX + aY). Then 
/ o o / _1 = is equivalent to saying that F o $ o and <& define the same 
morphism on P^ -> Pg. If F o $ o F _1 = (* (X F), *i(A, F)), then this means 

$o(A, F) - $!(A, F)* (X, F) = 0. (3.1) 

The expression on the left is a homogeneous polynomial of degree 2D in A and 
F whose coefficients are homogeneous polynomials in R[a,b, c, d]. So (|3.1[) gives 
2£> + 1 equations that cut out a closed subschemc of PGL2 defined over R. One 
checks readily that Aut^(S') is a subgroup of PGL2(<S') for every S. 

Next we argue that Aut^ is a finite group scheme over R when <f> has degree at 
least 2. The map Aut^ — > Speci? is quasi-finite. Indeed, it suffices to check this 
statement on geometric fibers, and Silverman has shown that Aut^(L) is a finite 
group for any algebraically closed field L [15] Prop. 4.65] Q 

Moreover, Aut^ is proper over Spec R. Indeed, since Aut^ and Spec R are noe- 
thcrian, this can be checked using the valuative criterion for properness using only 



Alternatively, £|4. II gives a slightly simpler proof of Silverman's result. 
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discrete valuation rings [HI Ex. II.4.11]. Let o be a discrete valuation ring with field 
of fractions k, and consider a commutative diagram 

Spec k s- Aut^ 



Spec o s- Spec R. 

The left vertical map is the canonical open immersion, and the right vertical map is 
the structure morphism. We must show there is a unique morphism Spec o — > Aut^ 
that makes the entire diagram commute. Without loss of generality, we may assume 
that R = o and that the lower horizontal arrow is the identity map. 

If v : k — > Z U {+00} is the canonical extension of the valuation on 0, then we 
may endow k with the structure of a non- Archimedean field by setting |x| = e~ v ^ 
for every x € k. (Note that we interpret e~°° as zero.) Since <p is defined over 0, it 
has good reduction. The Reduction Lemma asserts that every /c-automorphism of 
tfi also has good reduction. Equivalently, every k- valued point may be extended to 
an o-valucd point, which is what we wanted to show. 

We now know that Aut^ — > Spec R is a quasi-finite proper morphism. Zariski's 
main theorem tells us that it factors as an open immersion of i?-schemes Aut^ — > X 
followed by a finite morphism X — > Spec R. But Aut^ is proper, so any open 
immersion is actually an isomorphism. Hence Aut^ is finite over Spec R. This 
completes the proof of Theorem 11.11 



4. Algorithms 

4.1. Absolute Automorphism Groups — Method of Invariant Sets. Let 

F be an arbitrary field, and suppose <fi : P^ — > P^ is a morphism of degree at 
least 2. In this section we describe an algorithm to compute Aut^(F), where F is 
an algebraic closure of F. In fact, we will see that it also gives a field of definition 
E/F for the absolute automorphism group, although E is typically not the smallest 
such field. The idea is to use a finite Aut«^(i ;l )-invariant subset of P 1 (F) to produce 
a set of candidate automorphisms, and then to sort through the candidates by brute 
force. We begin by explaining Algorithm [TJ which determines the automorphism 
group if we already have such an invariant set T. Then we give a description of 
how one constructs T. 

Suppose that we are given a finite set T = {n, . . . , r„} C P J (F) with n > 3 on 
which Aut<f,(F) acts. Let E = F(ri, . . . ,r„) be the field of definition for T. Let 
s G Aut^E 1 ). Since we assume T is invariant, there must be a triple of distinct 
indices i,j, k 6 {1, . . . , n} such that s(ti) = Tj, sfa) = Tj, and s(tz) = r^.. 

Conversely, given a triple of distinct indices k <G {1, . . . there exists a 
unique element s e PGL2(-E) such that s(ti) = Tj, sfa) = Tj, and s(t3) = tu- One 
now determines if this candidate element s actually satisfies the functional equation 
s o tp — tp o s; if that is the case, then s £ Aut l j ) (E). See Algorithm [T] 

A natural candidate for an Aut<^(F)-invariant set is the set of fixed points of cj>. 
For let x G P 1 ^) be a fixed point, and let s G Aut^F). Then 



s(x) = s(cf>{x)) = 0(s(x)), 



(4.1) 
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Algorithm 1 — Compute Aut ( j ) (E) given an Aut^(i?)-invariant subset of V 1 (E) 
consisting of n > 3 points 

Input: 

• a nonconstant rational function <fi G E(z) 

• an Aut0(i?)-invariant subset T = {n, . . . , r„} C ¥ 1 (E), n > 3 
Output: the set Ant r j ) (E) 

create an empty list L 

for each triple of distinct integers i, j, k € {1, . . . , n}: 
compute s G PGL2(-E) by solving the linear system 

S(n) = Ti, s{T 2 )=Tj, s(t 3 ) = Tfc 

if s o (j) — (j) o s: append s to L 
return L 



which shows that s(x) is also a fixed point. A general rational function (f> of de- 
gree d > 2 has d + 1 > 3 distinct fixed points. Thus the choice T = Fix(0) will 
suffice in most circumstances. 

If instead Fix(</>) has cardinality 2, define T — cj)^ 1 (Fix(</>)). This set is Aut^(F)- 
invariant: if a; € Fix(0), y G _1 (a;), and s e Aut^F), then 

0(s(y)) = s(0(y)) = s(x). 

Since s(x) is a fixed point by equation (|4.1[) . we find s(y) e (Fix(</>)). Recall 
that we require #T > 3, and by construction we have #T > 2 since Fix(0) C T . 
If #T = 2, then T = Fix(^), and each of the fixed points is totally ramified for <f>. 
This implies that the derivative at each of the fixed points vanishes^ which in turn 
means each element of Fix(</>) has fixed point multiplicity 1. But the total number 
of fixed points of a map of degree d is d + 1 > 3, counting multiplicities, so we have 
a contradiction. (See, for example, [7J Appx. A].) 

Finally, suppose that Fix(</>) = {x}. We claim that > 2. For otherwise 

x is ramified for <j>, which implies that the derivative 4>'(x) vanishes there. But the 
fact that x is the unique fixed point of 4> means that in local coordinates centered 
at x our map is of the form z i— > z + a,d+iz d+1 + ■ ■ ■ with a^+i 7^ 0. So the derivative 
cannot vanish at x, and we must have #</) _1 (x) > 2 as desired. If 4t z 4>~ 1 {x) > 3, then 
set T = (f>~ l {x). Otherwise, 4r x {x) = {x,y}; set T = ^ 2 (x) = {x,y} U 4r x {y), 
which satisfies 3 < #T < d + 2. The argument in the preceding paragraph shows 
that T is Aut0(F)-invariant in all cases. 

Algorithm [5] gives the complete description of this method for computing the 
automorphism group of a rational function defined over an arbitrary field F. 

4.2. Method of Fixed Points. Let F be a finite field or a field of characteristic 
zero, and let <f> : — > a nonconstant morphism. For any (^-periodic point 
x G P 1 (-F), write per(a;) for its exact period — i.e., the minimum positive integer 
i such that 4> l (x) = x. If x is not periodic, write per(x) = +00. For each pair of 
integers i,j G {1,2}, define the following set: 

Z itj = {x G P^F) : per(a-) = i, [F(x) : F] = j}. (4.2) 



More precisely, the induced map T<f> on the tangent space TP^, is zero. 
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Algorithm 2 — Computation of Ant c / > (F) via the method of invariant sets 

Input: a field F and a rational function cf> G F(z) of degree at least 2 
Output: a field extension E/F and the set Aut^(i?) = Aut<p(F) 

if #Fix(0) > 3: 

compute Aut^E) using Algorithm [1] with T = Fix(0) and £ = F(T) 

else if #Fix(0) = 2: 

compute Aut</,(E) using Algorithm □ with T = (Fix(0)) and E = F(T) 

else: 

compute Aut ( J B) using Algorithm [1] with T = ^T 2 (Fix(0)) and E = F(T) 
return E and Aut r j,(E) 



We also define the following set of ordered pairs: 

W = {(*, y) : x G Zi,i, y G ^(x), : F] = 1}. (4.3) 

(More concretely, VF is the set of pairs of F-rational points such that x is fixed by 
<j) and (f>(y) = x.) These sets may be constructed by factoring the polynomials that 
define the fixed points of </>, the points of period 2, and the preimages of -F-rational 
points. We write for the set of unordered pairs of elements of a set Z. 

Suppose that s = ("5) A\it$(F) is nontrivial. The homogeneous polynomial 

defining the fixed points of s is jX 2 + (5 — a)XY — f3Y 2 . If s has a unique fixed 
point x, then F is a field of characteristic p > 0. (To see this, move the unique 
fixed point to infinity. Then s is a translation with finite order.) So we conclude 
that F = W q , that x is F-rational (because F is perfect), and that s has order p [5J 
Lem. 3.1]. Now 

s((/>(x)) = <fi(s(x)) = (f>(x), so that <f)(x) = x. 

Hence x G Choose u G PGL2(F) such that = 00; then usu -1 = (* \ ) 

for some A G F \ {0}. That is, s G it -1 ( 1 Fx / ^ ) u. In order to find all elements 
of Aut,p(F) of order p, it suffices to apply this technique to every x in the set Z\^\. 
See the first for-loop of Algorithm [3] 

Now suppose that s G Aiit<f,(F) has precisely two distinct fixed points x\ and x^. 
Then s(4>(xi)) = 4>(s(xi)) = 4>(x\), so that 4>{x\) G {aii,:^}- There are three 
possible cases: (1) (f> fixes both x\ and X2', (2) (/) swaps x\ and X2; or (3) 4>(x±) = X2 
and fixes X2 (perhaps after interchanging x\ and x^)- Since </> is defined over F, 
all Galois conjugates of a fixed point must also be fixed points. Thus in cases (1) 
and (2), either x\ and X2 are both F-rational, or they are quadratic conjugates over 
F. In case (3), both X\ and X2 must be F-rational. 

(2) 

If xi and X2 are both F-rational in case (1) — so that (xi, X2) G Z\ [ — then we 

may select u G PGL2(-F) such that u{x\) = 00 and ^(^2) = 0. Then usu^ 1 = ( 1 ) 

for some root of unity ( E F. If £ has order re, then re | d(d 2 — 1) by Proposition l2.2l 

Let T be the set of roots of unity in F that have order dividing d(d 2 — 1). We 

(2) 

loop over all distinct unordered pairs of elements (xi,X2) G Z\(, and check which 
elements of u _1 ( T 1 ) u lie in Aut^(F). See the second for-loop of Algorithm [3] In 
fact, this strategy works in case (2) when Xi,X2 are both F-rational, and in case 
(3). These correspond to looping over pairs (x\,X2) in Z\ ( and in W, respectively. 
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Now suppose x\ and X2 are quadratic Galois conjugates over F in case (1), so that 
(xi,x 2 ) G Then E = F(xi,x 2 ) is a quadratic extension. In this case, we may 

choose u G PGL^-E 1 ) satisfying u(x\) = oo and u(x2) = 0, so that usu^ 1 = ( i) 
for some root of unity £ G E. In particular, if the fixed points of s are a we 
may take 

A -a-bVd\ _ x (i \ (a+^^ -a* + b 2 d \ 

U ={l -a +b Vd)> S0S = U { l) U= [ I -a + ^j- 

Let a be the nontrivial clement of Gal(E/F). Since s G PGL 2 (-F), we have 

( (i + QbVd \ -(i + ^)bVd _ (1 + pbVd 

a { (i-o ; " (i-r) " (i-o ' 

which implies that = 1. If a acts trivially on £, we must have £ = — 1 (since s is 
nontrivial). Otherwise, E = F(£). In particular, F{x\,X2) is always a cyclotomic 
extension of F. 

If F = F q is a finite field, then £ G F g 2. Since the automorphism a acts by 
Frobenius on F g 2, it follows that = = 1. That is, a nontrivial element of 
PGL2(F 9 ) with two quadratic conjugate fixed points necessarily has order dividing 
q + 1. In this case, let A C F^ 2 be the unique subgroup of order q + 1; note that A 
always contains —1. 

If F is a field of characteristic zero, let C(X) = X d( - d2 -^ - 1. In this case, let A 
be the set of roots of the quadratic factors of C(X) over F together with — 1. 

(2) 

To detect s, loop over all Galois conjugate pairs in Z\ 2 and check which elements 
of u" 1 ( A 1 ) u lie in Aut^(F). See the third for-loop of Algorithm [3] The same 
strategy also applies if we are in case (2) and x\ and xi are quadratic conjugates. 

Remark 4.1. The first for-loop in Algorithm [3] will not terminate if F is an infinite 
field of characteristic p. But the remainder of the algorithm proceeds without 
modification and computes {s G Aut ( f ) (F) : gcd(order of s,char(F)) = 1}. 

Remark 4.2. When F is a number field, the technique in the proof of Proposition ll.4l 
gives further restrictions on the set of roots of unity one must include in the set A. 
It suffices to consider irrational and F-quadratic roots of X m — 1, where m is the 
product of all integers n such that <p(ri) \ 2[F : Q]. Here ip is Euler's function. 

4.3. Chinese Remainder Theorem Method for Number Fields. In this sec- 
tion, we assume that our field AT is a number field. As discussed in the previous 
section, Algorithm [3] computes Aut^AT). This algorithm requires computing the 
irreducible factors of degree at most 2 of a degree d 2 + 1 polynomial, namely the 
polynomial obtained by clearing denominators in the equation <fi o (f>(z) = z. As 
the degree becomes large, this quickly becomes impractical on computer algebra 
systems, even over K = Q. Thus, we have a need for an alternative algorithm over 
number fields. 

We use an approach that is ubiquitous in number theory: first compute the 
automorphism group over a residue field F„ for some finite place(s) v, and then use 
the local information to obtain a global answer. More precisely, our method is as 
follows. 

Let v = vq be a finite prime such that the reduction map Aut^(AT) — > Aut^(F^) 
is injective, e.g. v S as defined in Proposition [LU Now compute Aut^(F„) using 
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Algorithm 3 — Computation of Aut^(F) via the method of fixed points 

Lnput: a field F, finite or of characteristic zero, and 4> £ F(z) of degree > 2 
Output: the set Aut^(F) 

if F = F q is a finite field: 
create a list T = F* 

create a list A of £ e F* 2 \ {i} with = i 
else: 

let C{X) = X d ^ 2 -^ - 1 

create a list T of F-rational roots of C(X) 

create a list A of roots of -F-quadratic factors of C(X) and — i 

create a list L = [z] 

create the sets Zij , W defined in equations (|4.2j) and (|4.3j) 

for x £ Zi t i : 

choose u £ PGL2(F) such that u(x) = oo 
for A e T: 

set s{z) = u^ 1 (u(z) + A) 

if s o (j) = <p o s: append s to L 

(2) (2) 

for each pair (x, y) with x ^ y in Z\ { U Z\ [ U W: 

choose u £ PGL2(F) such that u(x) = oo and u(y) = 
for C € T n {1}: 

set s(z) = u~ 1 (£u(z)) 

if s o <f> = o s: append s to L 

(2) (2) 

for each pair of Galois conjugates (x, y) in f UZ^ : 

choose u £ PGL2(F(x, y)) such that u(x) = oo and = 
for £ e A: 

set s(z) = u _1 (£u(z)) 

if s o cj) = (j) o s: append s to L 

return L 



Algorithm [3J For each element in Aut0(F„), let / £ Aut(P] f ) be a lift of minimal 
height and check if / £ Aut^(K). If every element of Aut^(F^) lifts to an element 
of Aut c/,(K), then we are done. Otherwise, we will repeat this process, with some 
minor modifications; we explain these modifications now. 

Let vi be finite prime outside of S U {vq} and compute Aut < ^(F^ 1 ). Let G C 
Aut(Pjj if y Y\ v •) ^ e a su bset that surjects onto Aut^F^) for each i; this step is 
called the CRT (Chinese Remainer Theorem) step. For each element of G, choose 
a lift / <G Aut(P^) of minimal height. If/o(/> = </>o / then add / to the list Auts. 
If Auts surjects onto Aut^F^) for any i, then we are done. If not, then choose 
another prime ^ S U {«o, ■ • ■ 5 Vi} and repeat. 

In order to make this method into an algorithm, we need to provide a terminating 
condition. Write N(v) for the norm of a finite prime v. We claim that if ]X N(uj) > 
(2M)<- K: ®' , for some explicitly computable constant M, then Auts = Aut^(if), even 
if Auts does not surject onto Aut^,(F Vi ) for any i. We will spend the rest of the 
section proving this claim via the theory of heights. 
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Let H k '■ P^A) — > K>i denote the relative multiplicative height for A and let 
L 2 (f) denote the L2-norm of a polynomial /. See, for example, [9J B.2, B.7] for 
definitions. 

Proposition 4.3. Let T,T' C P X (A') be Galois invariant sets of order at least 3, 
and let /t,/t' G K[ w i z ](o) be square-free polynomials such that V(fx) = T and 
V{f T >) = T'. Then for any s G Aut(P^) C P 3 (A) such that s(T) = T' , we have 
H K (s)<ti KA ®L 2 (f T ) 3 L 2 (fT>) 3 . 

Proof. Let s be as in the statement of the Proposition. Let Ti,T2, T3 be 3 distinct 
elements of T, and let rji := s(rj) G T' . In coordinates, we write n = (t^o : Ti,i) 
and rji = (77^0 : Vi.i)- Since an automorphism of P 1 is determined by its action on 
3 elements, we have an expression for s = ("5) m terms of Ti t j,r]i t j, i.e. 

a = E CTe S 3 (sgll cr)S . (1) C CT ( 2 )D .(3), P = E CTG S 3 ( S S n cr ) j4 <r(l) C ' CT (2)£' ( T(3), 
7 = E CTe s 3 ( s g n ct )4t(i)£<t(2)At(3), <5 = E CTe s 3 ( s g n cOAxU) ^(2)^(3), 

where ^ = rj,o»7i,i, 5* = -Ti^i, Cj = -r^.o^.O) an d A = ^,1^,0- 

This expression allows us to obtain a bound on the local height of s. Let i> be 
any place of A and let e v = 6 if v | 00 and e — 1 if v \ 00. Then, by the triangle 
inequality, 

\ot\v < Sv ■ max |5 cr (i)C , (T (2)L» cr (3)| 1 , < [ max{|r i0 | w , |rii|t>} • max{|??iok l^ilf}- 

l<i<3 

One can easily check that the same bound holds for |/3|,j,|7|„,|(>|.„. It follows that 
Hk{s) = JJmax{|a|„, \p\ v , \y\ v , \6\ V }^ K ^ 

V 

< l[e^ ■ [] niaxiKoL \m\ v }^^ ■ max^U \va\ v } M] 

v l<i<3 

= 6 [^:Q] Y[ H K (Ti)H K (m). 
l<i<3 

Since Hk(ti) < ^(/t) and HxiVi) < L 2 (fx') [SI Lemma B.7. 3.1], this completes 
the proof. □ 

Corollary 4.4. Lei <f> G A(z) &e a rational function of degree > 1, let T C P 1 (A) 
be the Galois invariant set constructed in Algorithm® {with F = A), and let fx be 
a square-free polynomial such that V(fx) = T. Then every element 0/ Aut^(A) C 
P 3 (A) has relative multiplicative height bounded by 6^ K:II ^L 2 (fx) e 

We will take this height bound 6^ K:<! ^L 2 (fx) ( ' to be our explicit constant M. Now 
we need to show that if ni N ( v > {2M) 2 ^ K '- Q \ then every element of Aut^(if) is 
a lift of an element of Aut^(F I)i ) of minimal height. We will need the following 
two lemmas. 

Lemma 4.5. Let b be a nonzero fractional ideal of Ok, and write it as a quotient 
b = b + /b~ of relatively prime integral ideals. Then Hk(/3) > N(b + ) for all nonzero 
/3 Gb. 
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Proof. Since (3 £ Ok, we have |/5 
HkW) = J] max{l, 

v\oo 



\ v < 1 for any finite place v. Therefore 



lv] [K v :Q v] " max{1) ^ 



[A„ 



v\oo 



[A\, 



n i 

?)(b)<0 



u(b)<0 
|[it„:Q„] . 



n i 

ujoo 
i>(b)>0 



where the last equality follows from the product formula. Let e p be such that 



b = J]p e p. Since /? € b, v(f3) > e Pv , so 



-[iC*:< 



vS > N( Pt ,) e ^ 



□ 



Lemma 4.6. Let a C Ok be an integral ideal, and let p a : ¥ 71 (Ok) — > ^(Ok/o) 
denote the canonical projection. For each j3 = (/3q : /?i : ■ • • : f3 n ) <E P"(Oif/a), 
£/iere is at most one element a = (ao : tt\ : ■ ■ ■ : a n ) 6 P^iP) with Hk(o) < 

(2-[^QlN(a)) 1/2 . 

Proof. Let a, a' € P n (0 K ) be such that H K {a),H K {a') < (2^ K: ^N(a)) 1/2 and 
such that p a (aO = Pa^')- Since a € P"(Ck), there exists a coordinate io such that 
cti a a. It follows that oi i $ a too. 

Then for each i and each place v, an argument as in Proposition 14.31 shows that 



1. 



Oijr. 



Hi 

a' 



< 2 max 



Of 



max ■ 



Taking the product over all v gives H K (-^- - < 2^ K ^H K {a)H K (a'). The 

latter is less than N(o) by hypothesis, and — ^- lies in the fractional ideal 

(tti a^ so the preceding lemma implies that — = That is, a — a' . □ 

Proposition 4.7. Let uqi • • • j Vn be finite places of K such that 

(1) the reduction map Aut^(A") — > Aut^F.^) is injective for all i, and 

(2) Il l N (^) > 2l if:Q ]M 2 , where M = Q^ K ^L 2 {f T f is as in Corollary^ 
For any tuple (gi) € Yii Aut^F^), let gx G Aut(P^) be a simultaneous lift of each 
gi of minimal height. If ga ^ Aut^AT), then (gi) ^ im(Aut0(A) — > J\ i Aut0(F„J) 

Proof. Assume that (g^ e im(Aut0(A) FJi Aut^F.^)) and let g' e Aut^(A) 
denote its pre-image. (The automorphism g' is unique by assumption (1).) By 

CorollaryOl H K {g') < M < (2^ K: ^ IL N (>i)) 1/2 - B Y LemmaSH g 1 must have 
minimal height among all lifts, so g' = gx G Aut^(iT) □ 

There are a few technical details that we have left out in our description of 
Algorithm 21 mostly in the step where we decide whether to terminate and in 
the Chinese Remainder Theorem step. These details allow us to avoid extraneous 
computation. We give an example here and the curious reader can find the rest in 
the source code. 

It is possible for the reduction of Aut^(A) to be a proper subgroup of Aut^(F„) 
for all places v of good reduction. Consider the rational map 4> = 2z 5 . One can use 
the method of invariant sets to check that 



Aut rf 



[z,iz 



-iz, (V2z)-\i{V2z)-\ ~(V2z)-\ -i^z)- 1 } , 
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Algorithm 4 — Computation of Aut^-ftT) via the Chinese Remainder Theorem 

Input: a number field K and a rational function <f> G K(z) of degree d > 1 
Output: the set Aut^A') 

choose T as in Algorithm [2] and set M = 6^ K ^L 2 (f T ) 6 
create an empty list L, and set a = (1) 

for v a prime of good reduction at v such that Aut<p(K) — > Aut0(F„) is injective: 
compute Aut^(F„) using Algorithm [3] 
if Aut (F„) = {z}: 

return {z} 
else: 

append Aut^(F w ) to L, and set a = ap„ 
Set L' = CRT(L) and initialize an empty list Auts 
for ,s in L': 

set s' G PGL2(Ok) to be a lift of s of minimal height 
if H K (s') < M and s' o = o s': 
append s' to Auts 
if N(a) > 2^ K: ^M 2 or if #Auts = # Aut (F,„) for any v \ a: 
return Auts 



which is a dihedral group of order 8. For all primes p > 2, at least one of —1, 2, —2 
is a square in F p . Therefore, Aut^,(F p ) always contains Z/2 x Z/2 or Z/4 as a 
subgroup. As the algorithm is stated, we would compute a lift of every element in 
11^5 Aut^(Fp). However, by p — 7 one can already recognize that Aut^(Q) C Z/2 
since Aut0(F 5 ) = Z/4 and Aut^(F7) = Z/2 x Z/2. Our code checks for group- 
theoretic properties like this when deciding whether to terminate. 

It is important to build in as many early termination conditions as possible, 
since typically the elements of Aut^(K) have significantly smaller height than the 
theoretical bound. This is of course true when Aut^(-R') is trivial, but it remains 
true even in the non-trivial case. For example, consider the functions in the last 3 
lines of Table HJ The height bound for </> = 345025251z 6 is over 50 digits, and the 
height bounds for the other two are over 100 digits. In contrast, the heights of the 
automorphisms are (from last to first) 2601, 101 and 11. 

5. Examples 

In this section, we compute some examples to give an idea of the running time. It 
is hard to produce "random" examples that have non-trivial automorphism group. 
Therefore, we present some hand-selected examples with non-trivial automorphism 
group which demonstrate the correctness of the algorithm (see Table [1) . Then 
we present median running times for randomly generated rational maps of varying 
degrees and varying heights (see Table [2]). All of the randomly generated functions 
had trivial automorphism group. 

Our computations indicate that the fixed-point method is faster for rational 
functions of small degree, but that the CRT method is a better choice once the 
degree is larger than 15. In our implementation, the main bottleneck in the fixed 
point algorithm is in computing Z\,i and Z2.2; this requires computing the quadratic 
factors of a degree d 2 + 1 polynomial . It is possible that the fixed-point method 
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<p 


CPU 1 
CRT 


ime 

FP 


Aut^(Q) 


2 

z — 22 — 2 
Q-,2 2^- l 


0.338825 


0.058 


/ \ ±l 
z ' i+T 


2 2 -42-3 


0.101023 


0.020238 


-2-1 -1 


— 3z' 2 — 22+2 


' 2 5 Z+l 


2 — 52 +102 — 52 

— 52 4 + 102 ii — 52 + 1 


2.445185 


0.205686 


- £ - 1 1 1 22-1 ~ Z + 2 
7 2—1' 5 2 7 2 — 2 ' 2+1 5 

2 — 2 2 + 1 —1 2—1 —2 — 1 22—1 
22 — 1 7 22— I 7 2 — 1 1 2 1 2 — 2 7 2 + 1 


2 5 -202 4 +302 3 + 102 2 -202+3 


0.583004 


0.035206 


„ 2-2 -1 2-1 -2-1 22-1 


— 02 — 02 +4U2 J — OU2 — 02+4 


7 Z2 — 1 2 — 1 ' 2 ' Z — 2, ' Z-\-L 


32 2 -l 

2 y — 32 


0.227719 


0.042237 


2' \ Z-\-l J \ 2 — 1 / 


2 3 -2l2 2 -32+7 


0.648 


0.024699 


„, —1 2—1 —2—1 


-72 i! -32 2 +2l2 + l 


Z 1 2 ' 2 + 1 ' 2-1 


2 11 +662 6 -ll2 
-ll2 1 "-662 5 + l 


0.676026 


0.069261 


2,-1/ Z 


32z 5 


0.681829 


0.035364 


±z 


4> = (z-™y f = =^ 


1.663506 


2.962213 


„ -32+5 




112 + 3 


0=(z 19 )/;/ = ^2_2 


43.554845 


4.678652 


„ -172-7 -3z+5 -132-53 
> -52+17' 112+3 ' 1012+13 


345025251z 6 


337.721855 


0.073642 


z,l/(260U) 



Table 1. Running times for automorphism_group_QQ on rational 



functions with non-trivial automorphism group. 



can be made feasible for larger degrees by implementing a faster method for finding 
quadratic factors of large degree polynomials. 

These examples were computed on a Macbook Air (Apple, Inc.) running Mac 
OS X 10.7.2 with a 2.13 GHz Intel Core 2 Duo processor and 2GB of RAM. They 
were run with Sage 4.7.2 which was released on October 29, 2011. All running 
times are listed in seconds. 

6. Conjugate Rational Functions 

Let F be a number field or finite field, and let <p, ip : Pp — > Pp be a pair of 
endomorphisms of the same degree d. We return now to the question presented in 
the introduction: 

Does there exist a change of coordinate / £ PGL^i* 1 ) such that 

ip = f o <p o / ; i.e., are cj> and ip conjugate over Fl 
This question can be dealt with both theoretically and computationally in a manner 
similar to that of automorphism groups. We briefly sketch the theoretical side; the 
proofs are straightforward modifications of arguments presented earlier. 

Definition. Fix a nonnegativc integer d, and let </>, ip : — > P^ be two endomor- 
phisms of degree d. The conjugation scheme of the pair (<p, ip) is the i?-scheme 
Conj^, ^ represented by the functor Conj^ : i?-Alg -4 Grp defined by 

Conj . JS) = {.f £ Aut(P^) : / o o = V+ 
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d 


50 


10 2 


Height Bound 
10 3 10 4 


10 5 


10 6 


Fixed Point 


3 
6 
9 

12 
15 


0.0124450 
0.0254290 
0.0773140 
0.282548 
0.980427 


0.0133855 
0.0263755 
0.0808070 
0.2629635 
1.0104325 


0.0133215 
0.0270010 
0.0807205 
0.2917585 
1.0165010 


0.0132855 
0.0270970 
0.0908770 
0.2940025 
1.0359760 


0.0127770 
0.0294735 
0.0938780 
0.3132295 
1.1184880 


0.0129115 
0.0307725 
0.0961555 
0.3472265 
1.1377990 


CRT 


10 
15 
20 
30 
40 


0.1932860 
0.5662205 
1.8513510 
17.702872 
106.19608 


0.1968665 
0.5173635 
1.9718745 
18.402277 
102.52844 


0.2083450 
0.5148960 
1.9562015 
17.750321 
106.19090 


0.2070960 
0.5360245 
1.7734270 
20.359524 
113.88063 


0.2180315 
0.4533855 
1.9558800 
18.483440 
113.06002 


0.2347400 
0.5041420 
1.9808705 
18.752402 
122.07266 



Table 2. Median running times for the fixed point and CRT al- 
gorithms on 100 random rational functions with given degree and 
height bound. 



Theorem 6.1. Let R be a commutative ring, let d > be an integer, and let 4>,ip : 
Pjj — > Pjj be endomorphisms of degree d. Then the functor Conj^ is represented 
by a closed R-subscheme Conj^^ C PGL2. If moreover d > 2, then Conj^^, is 
finite over Spec R. 

The theorem does not preclude the possibility that Conj^ ^ is the empty scheme; 
in fact, it is typically empty when d > 2. The group scheme PGL2 has relative 
dimension 3 over R, and the space Rat^ of endomorphisms of P 1 of degree d has 
relative dimension 2d + 1 > 3 over R. So for a fixed </> € Rat £ ;(i?), a general choice 
of ip will yield Conj . ^ = 0. 

Remark 6.2. When Conj^, ^ is not the empty scheme, it is a principal homogeneous 
space for Aut^ (or Aut,/,). 

In the case d > 2 of the theorem, in order to establish that Conj^ ^ is finite over 
Spec R, one must argue that it is proper and quasi-finite. Properness follows from 
a direct generalization of the Reduction Lemma. Quasi-finitcness may be proved 
by taking R = F to be an algebraically closed field and using a variation of our 
method of invariant sets in £14.11 The basic idea is to replace the invariant set T 
with two sets T^, C F 1 (F) such that 

• > 3, and 

• s(T ) C T. 4 , for every s E Conj^(F). 

The sets T^,Ty, may be taken as the fixed points of <p and ip, respectively; if there 
are not enough fixed points, then one may use pre- images exactly as in Algorithm!]] 
We leave the details to the reader. 

Finally, we note that the algorithms in H4.ll and £14.31 may be adapted (as in the 
previous paragraph) to give algorithms for computing Conj^ ^(F) and Conj^, AF) 
when F is a finite field or number field. Again, we leave the details to the interested 
reader; note that the main technical tool — Proposition 14.31 — applies in this 
more general situation. We close with the remark that, in particular, this strategy 
gives an algorithmic procedure for determining if Conj^, ^ is the empty scheme or 
not. More concretely, this means that one can determine effectively if two rational 
functions <f> and tp are conjugate. 
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